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Abstract 

In  this  thesis,  a  class  of  time  series  models  for  forecasting  a  hurricane’s  future 
position  based  on  its  previous  positions  and  a  generalised  model  of  hurricane  motion 
are  examined  and  extended. 

Results  of  a  literature  review  suggest  that  meteorological  models  continue  to 
increase  in  complexity  while  few  statistical  approaches,  such  as  linear  regression, 
have  been  successfully  applied.  An  exception  is  provided  by  a  certain  class  of  time 
series  models  that  appear  to  forecast  storms  almost  as  well  as  current  meteorological 
models  without  their  tremendous  complexity. 

A  suggestion  for  enhancing  the  performance  of  these  time  series  models  is 
pursued  through  an  examination  of  the  forecast  errors  produced  when  these  models 
are  applied  to  historical  storm  tracks.  The  results  uncover  no  patterns  that  can  be 
exploited  in  developing  an  improved  model  and  suggest  that  there  are  meteorological 
processes  or  factors  at  work  beyond  those  that  can  be  modeled  with  the  available 
historical  data  base. 

The  statistical  structure  of  the  time  series  approach  is  exploited  to  develop  a 
practical  method  for  determining  prediction  regions  which  probabilistically  describe 
a  hurricane’s  likely  future  position.  The  Monte  Carlo  approach  used  to  develop  these 
prediction  eUipses  is  seen  to  be  applicable  for  predicting  areas  subject  to  risk  from 
hurricane  landfall. 
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Developing  Prediction  Regions  for  a  Time  Series 
Model  for  Hurricane  Forecasting 


I.  Introduction 

Hurricanes  Andrew  (1992)  and  Hugo  (1989)  are  well  known  recent  examples 
of  the  tremendous  destructive  capability  of  hurricanes.  The  difficulty  that  forecast¬ 
ers  had  predicting  the  time  and  location  of  landiall  is  remembered  less,  but  had 
significant  ramifications  for  those  areas  prepared  and  unprepared  for  the  storms. 
Hurricane  Emily  (1993)  is  not  remembered  as  well  because  it  did  not  make  landfall, 
but  forecasters  were  uncertain  of  this  until  it  curved  eastward  and  moved  away  from 
land.  These  examples  illustrate  the  need  for  accurate  and  timely  hurricane  forecasts, 
especially  near  the  coast. 

Hurricanes  are  part  of  a  class  of  storms  generally  called  tropical  cyclones  and 
tropical  storms.  In  the  Atlantic,  they  form  off  the  west  coast  of  Africa  in  warm 
equatorial  waters  from  about  June  1  to  November  1.  The  first  stage  of  development 
is  called  a  tropical  depression,  when  an  area  of  organized  circulation  forms.  Tropical 
depressions  are  characterized  by  poor  boimdary  definition  and  sustained  wind  speeds 
of  less  than  34  knots.  The  next  stage  is  tropical  storm  and  is  marked  by  an  organized 
circulation  pattern.  Tropical  storms  are  characterized  by  sustained  wind  speeds 
from  34  to  63  knots.  The  last  stage  is  hurricane.  A  well  defined  three  dimensiomd 
circulation  structure  is  evident  and  sustained  wind  speeds  are  64  knots  or  greater. 
Diminishing  storms  pass  through  these  stages  in  reverse,  usually  after  landfall  or 
after  they  reach  northern  ocean  waters  too  cool  to  add  enough  energy  to  sustain 
them. 
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Hurricanes  are  not  differentiated  farther  in  this  thesis.  The  te.'ms  hurricane, 
storm,  and  tropical  itorm  will  be  used  to  mean  all  storm  stages  generally,  when  not 
otherwise  obvious  in  context. 

The  National  Hurricane  Center  (NHC)  is  the  agency  tasked  with  making  official 
tropical  cyclone  forecasts.  During  an  active  storm,  they  provide  six  hour  update  of 
forecasts  out  to  72  hours  (twelve  six-hour  periods).  Hurricanes  have  been  hard  to 
forecast  because  they  do  not  consistently  follow  the  same  patterns  that  other  weather 
systems  do.  The  motivation  for  wanting  to  predict  their  motion  is  to  minimise  the 
damage  from  unexpected  landfall. 

A  certain  amount  of  position  and  intensity  uncertainty  arises  from  the  scarcity 
of  data  about  a  particular  storm,  especially  early  in  its  formation  and  in  remote  areas 
of  the  ocean.  "Best  track”  data  is  usually  determined  after  the  storm  has  ended  and 
represents  the  best  estimate  of  storm  position  developed  from  all  available  sources, 
such  as  satellite  imagery  and  aircraft,  ship,  and  soimding  reports.  Best  track  data 
is  used  when  developing  forecast  models  and  testing  their  performance. 

Over  the  past  few  years,  the  number  of  aircraft  ffights  to  get  real-time  position 
and  other  data  has  steadily  decreased,  primarily  because  of  reduced  budgets,  not  a 
redaction  in  the  number  of  storms.  Forecasters  rely  more  and  more  on  satellite  im¬ 
agery  for  real-time  information  about  a  hurricane,  including  position,  intensity,  and 
wind  speed.  Appendix  A  includes  documentation  of  current  methods  to  determine 
position  and  intensity  from  satellite  data. 

In  this  thesis,  attention  is  focused  on  a  class  of  time  series  models,  developed 
by  Curry  (3)  and  extended  by  Mott  (12),  that  forecast  hurricane  position  based 
on  a  generalized  model  of  hurricane  motion  and  a  particular  hurricane’s  previous 
positions.  The  specific  objectives  of  this  research  are  to  exploit  this  researcher’s 
meteorological  background  to: 
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(i)  review  the  cttnent  literature  to  summarise  the  current  "state  of  the 
art”  in  hurricane  forecasting; 

(tt)  develop,  based  on  published  results,  an  assessment  of  the  relative 
performance  of  current  methods  for  hurricane  forecasting; 

(in)  extend  Mott’s  and  Curry’s  efforts  by  developing  a  practical  method 
for  determining  "prediction  regions”  which  describe  a  hurricane’s  likely 
future  position;  and 

(ii».<  assess  the  viability  of  improving  Curry’s  and  Mott’s  approaches  by 
^YRtnining  those  storms  that  appear  in  some  way  "anomalous”,  that 
is,  not  well  forecast  by  their  models. 

This  thesis  first  examines  the  work  accomplished  by  Curry  and  continued  by 
Mott.  Pertinent  aspects  of  the  model  are  discussed  followed  by  a  literature  review 
covering  other  statistical  and  meteorological  models  developed  since  Curry’s  review 
in  1985.  The  models  are  then  compared  on  the  basis  of  their  reported  results. 

An  important  contribution  of  this  thesis  is  the  adaptation  of  some  theoreti* 
cid  results  developed  by  Curry  that  enable  one  to  develop  prediction  regions  for  a 
hurricane’s  future  position  which  will  help  quantify  the  uncertainty  associated  with 
a  forecast  (3).  In  particular,  a  Monte  Carlo  approach  is  developed  which  sidesteps 
many  of  the  complexities  inherent  in  Curry’s  proposed  methodology  and  which  ad¬ 
ditionally  has  the  potential  to  be  used  to  develop  landfall  prediction  regions. 

This  is  followed  by  an  analysis  of  the  forecast  errors  when  Mott’s  model  is 
applied  to  historical  data.  In  particular,  those  storms  that  appear  to  be  potential 
outliers,  producing  much  larger  than  average  forecast  errors,  are  identified  and  ana¬ 
lyzed  to  see  if  they  share  any  common  characteristics  (such  as  being  early  or  late  in 
the  season)  that  could  be  easily  incorporated  into  the  time  series  models  developed 
by  Curry  and  Mott. 
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II.  Time  Series  and  Hurricane  Forecasting 

Box-Jaikms  time  series  models  provide  one  way  to  forecast  random  variables. 
They  are  distinguished  from  other  modek  in  their  use  of  a  combination  of  past  values 
and  observational  errors  to  explain  a  present  value  or  make  a  forecast.  First  Curry, 
then  Mott  used  this  approach  to  forecast  hurricane  tracks  with  &vorable  results.  We 
review  these  models  in  the  first  three  sections  of  this  chapter. 

Other  statistical  approaches,  such  as  linear  regression,  have  been  applied  to 
hurricane  forecasting,  generally  with  less  success.  Along  with  statistical  models,  the 
NHC  uses  a  number  of  meteorological  models  to  forecast  hurricanes  in  an  attempt  to 
minimize  the  error  from  any  one  model  (17).  Some  of  these  recent  models,  including 
NHC90,  are  discussed  in  later  sections  of  this  chapter.  The  chapter  concludes  with 
a  comparison  of  the  forecast  errors  published  by  Curry,  Mott  and  the  NHC. 

2.1  Box- Jenkins  Time  Series  Models 

The  fundamental  basis  for  both  Curry’s  and  Mott’s  time  series  models  of  hurri¬ 
cane  motion  are  the  autoregressive-moving  average  (ARMA)  models  developed  and 
explained  by  Box  and  Jenkins  (1).  In  an  ARMA  model,  a  random  variable,  1^,  is 
expressed  as  a  weighted  series  of  its  values  at  past  times,  a  weighted  series  of  past 
errors,  and  the  random  “white  noise”  error  associated  with  the  current  observation. 
The  white  noise  errors  are  assumed  to  be  independent  and  normally  distributed 
with  mean  of  zero  and  constant  variance.  An  ARMA  model  is  built  empirically,  but 
knowledge  of  the  generating  process  sometimes  helps  in  selecting  the  best  model. 
The  following  sections  describe  the  various  types  of  univariate  ARMA  models  where 
a  random  variable  is  explained  only  by  its  past  values  and  its  associated  rauidom 
errors. 
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S.1.1  Autoregressive  Series.  An  autoregressive  series  has  the  form: 


Yt  =  ^Yt~\  +  <i>iYt-i  + - 1-  ^Yt-p  +  ej  +  C  (1) 

where  the  iff  coeffidents  are  the  relative  weights  apphed  to  past  observations,  et  is  a 
random  error,  and  C  is  a  constant  for  the  process.  The  subscript  t  u  the  time  index. 
The  number  of  past  observations,  p,  used  in  equation  (1)  is  referred  to  as  the  order 
of  the  AR  series. 

2.1.2  Moving  Average  Series.  A  moving  average  series  has  the  form: 

Yt  =  et  —  6iet~i - 6t-q^t-q  +  C  (2) 

where  the  e’s  are  the  random  errors  that  occurred  at  previous  times,  the  0 ’s  are 
their  relative  wdghts,  and  C  is  a  constant  for  the  process.  The  MA  is  of  order  q  and 
expresses  Yt  as  a  function  only  of  current  and  previous  errors. 

2.1.S  ARM  A  Models.  An  ARMA  model  incorporates  elements  of  both  AR 

and  MA  series  and  has  the  form: 

Yt  =  if>iYt-i  +  if>2Yt-2  +  •  — <f>pYt-p  +  e*  —  6iet-i  —  —  •  •  •  —  OgCt-q  +  C  (3) 


and  is  said  to  be  of  order  (p,q). 

The  order  of  an  ARMA  process  is  diffic’ilt  to  determine  solely  from  a  series 
of  observations.  Some  knowledge  of  the  underlying  process  is  usually  helpful.  Most 
ARMA  models  can  be  adequately  described  by  order  (2,2)  or  less,  even  those  which 
are  generated  by  higher  order  processes  (9:362).  Examining  the  autocorrelation  and 
partial  autocorrelation  functions  of  a  postulated  model  will  usually  help  indicate  the 
order  of  2m  ARMA  series.  The  autocorrelations  give  a  quantitative  indication  of  the 
interdependence  between  observations  in  the  series  while  the  partial  autocorrelations 
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provide  mfoimation  about  the  order  of  the  regressive  component.  Fitting  an  ARM  A 
model  to  data  is  traditionally  a  somewhat  complicated  process  which  requires  (i) 
identifying  the  order  of  the  model,  (ii)  estimating  its  parameters,  and  (iii)  assessing 
the  adequacy  of  the  fitted  model.  Makridakis,  et  al,  fpve  a  provide  a  good  description 
and  illustration  of  this  process  (9).  Several  software  packages  are  available  that 
estimate  the  parameters  of  time  series  models  from  a  series  of  data.  A  final  check 
of  the  adequacy  of  the  fitted  model  is  provided  by  determining  if  the  residual  errors 
are  statistically  the  same  as  white  noise. 

ARM  A  models  can  be  developed  only  for  series  that  have  a  constant,  or  sta¬ 
tionary,  mean  and  variance.  A  relatively  straightforward  way  to  induce  stationaiity 
in  a  non-stationary  data  series  is  to  “difference”  the  data  by  subtracting  each  ob¬ 
servation  from  its  successor.  The  series  of  differences,  or  lags,  is  then  used  in  the 
model  rather  than  the  raw  data.  Specifically,  lag  1  is  defined  as  observation  1  minus 
observation  2,  lag  2  is  observation  2  minus  observation  3,  and  so  on.  Sometimes 
more  than  one  difference  is  needed  to  induce  stationarity,  although  each  difference 
reduces  the  series  length  (and  degrees  of  freedom)  by  one  (10:45).  Box  and  Jenkins 
refer  to  this  aspect  of  time  series  modeling  as  integration,  with  the  resulting  model 
linearly  approximating  the  differenced  series.  Curry  found  that  the  first  differences 
of  hurriceine  position  were  sufficiently  stationary  to  model  with  ARMA  models. 

Z.  1.4  Transfer  Function  Models.  An  extension  of  univariate  ARMA  models 
is  to  include  past  values  of  one  or  more  other  explanatory  variables  in  the  model 
equation.  The  multivariate  case  has  the  form: 

Yt  =  +  +  (4) 

-f  UoXt  —  UiXt-i  —  •  •  •  —  UmXt-m 
+  CoZt  —  —  •  •  •  —  CnZt-n  +  Ct 
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where  X  and  Z  denote  the  other  random  variables,  u  and  ^  denote  their  respective 
weights  and  m  and  n  are  their  respective  orders. 

2.2  Curry ’3  Model 

The  purpose  of  Curry’s  research  was  to  develop  a  quick  and  accurate  means  of 
forecasting  hurricanes  so  that  fewer  people  would  have  to  evacuate  and  those  who 
did  would  have  greater  forewarning  in  the  event  of  hurricane  landfall.  Based  on 
their  relative  simplicity,  he  chose  a  time  series  approach,  eschewing  the  vastly  more 
complex  modeling  schemes  developed  and  used  by  the  NHC  (3:3). 

Gurry  proposed  a  first-difference,  bivariate,  fifth-order  AR  model  to  model 
hurricane  tracks.  In  his  model,  he  used  the  series  of  first  differences  of  latitude  and 
longitude  as  predictors  of  each  other.  He  further  refined  it  as  a  "threshold  model”  by 
dividing  it  into  7  latitude  bands,  or  thresholds.  The  model  has  distinct  coefficients 
in  each  band.  Forecasts  are  made  for  successive  6-hour  periods,  recursively  using 
each  forecast  as  a  current  position  through  72  hours.  In  each  latitude  band,  Curry’s 
model  has  the  form: 

LAt  —  LAt-i  =  —  LAt^2)  +  <l>ii,4{I'At-A  —  LAts)  (5) 

+  —  LOt-2)  +  —  LOt-s) 

+  Cl  +  au 

LOt  —  LOt-i  =  ~  +  <h\,\{LAt-\  —  LAt^^)  (6) 

+  4^1,\{,L0t-l  —  LOt-2)  +  <f>22,4{LOt-4  —  LOt-s) 

+  Cj  +  ajt 

where  LAt  and  LOt  denote  a  storm’s  latitude  and  longitude  at  time  t  and  the  a  terms 
are  random  errors.  Note  that  Equations  (5)  and  (6)  are  essentially  AR  transfer 
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function  models  with  the  same  order  applied  to  the  first  difference  series  of  both 
latitude  and  longitude. 

To  develop  this  model,  Curry  started  with  the  best  track  data  firom  storms 
1886-1983.  Best  track  data  is  determined  after  all  data  sources  report,  usually  after 
a  storm  is  finished,  and  represents  the  best  estimate  of  the  actual  path  of  a  storm. 
For  accuracy,  he  included  only  portions  of  storms  of  tropical  storm  strength  (34  knot 
wind  speed)  or  greater,  storms  in  the  Atlantic  basin,  and  those  that  occurred  after 
1945.  Position  reports  following  landfall  on  the  continental  U.S.  were  also  deleted 
since  a  storm’s  behavior  over  land  can  be  expected  to  be  different  from  that  over 
water  (3:134).  The  final  data  base  contained  362  storm  tracks  vdth  approximately 
10,000  position  reports. 

Curry  differenced  this  data  after  he  assumed  and  confirmed  that  the  storm 
tracks  were  non-stationary.  That  is,  hurricane  tracks  clearly  do  not  have  a  constant 
positional  mean;  they  move  across  the  ocean  during  their  life  and  their  positions  do 
not  vary  randomly  about  a  mean  latitude  or  lon^tude.  However,  the  first  difference 
series  of  storm  position  appears  to  reasonably  fit  the  criteria  of  mean  and  variance 
stationaiity  for  time  series  modeling.  That  is,  the  change,  or  lag,  in  position  from  one 
period  to  the  next  appears  to  hover  around  a  constant  mean  value.  Curry  identified 
nine  stationaiity  categories  based  on  the  behavior  of  the  first  and  second  differences 
of  latitude  and  longitude  (3:37).  He  worked  only  wi^h  those  storms  that  appeared 
to  be  “velocity”  stationary  under  the  assumption  that  storms  in  the  other  categories 
were  not  sufficiently  stationary  to  be  modeled  with  time  series  methods. 

After  further  examining  the  autocorrelation  and  partial  autocorrelation  func¬ 
tions  of  successively  higher  order  AR  models,  Curry  found  that  an  autoregressive 
transfer  function  model  of  order  5  best  explained  a  storm’s  present  position.  The 
second,  third,  and  fifth  lags  were  generally  weak  and  so  were  dropped  from  his  final 
model  (3:39),  as  seen  in  Equations  (5)  and  (6). 
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To  obtaiB  forecasts  firom  this  model,  a  one-step  ahead  forecast  position  is  based 
on  the  five  most  current  position  reports.  To  obtain  a  two-step  ahead  forecast,  the 
one-step  ahead  forecast  is  treated  as  the  last  ohaerved  position  and  is  used  along  with 
the  four  most  recent  position  reports  to  obtain  a  pseudo  one-step  ahead  forecast. 
Forecasts  for  lead  times  up  to  twelve  steps  ahead  are  computed  in  a  similar  manner 
(3:65). 

Since  the  general  motion  of  a  storm  changes  with  latitude,  Curry  subdivided 
the  model  into  latitude  bands  in  which  the  model  coefficients  (^’s)  differ.  When  a 
forecast  crosses  a  latitude  threshold,  coefficients  for  the  new  band  are  used  for  all 
forecasts  in  the  new  band.  Great  circle  distance  between  the  forecast  and  best  track 
positions  is  reported  as  the  forecast  error  in  nautical  miles  (nmi),  or  kilometers  (km) 
(17:526).  The  same  error  is  reported  by  the  NHC.  Curry’s  model  and  coefficients 
are  included  in  Appendix  B. 

To  test  the  validity  of  his  model,  Curry  calculated  the  model  coefficients,  then 
singly  removed  each  storm  from  the  data  base  and  recalculated  the  model  coefficients. 
Hypothesis  tests  showed  no  statistical  differences  in  the  coefficients  were  introduced 
by  any  one  storm.  Finally,  he  validated  the  model  with  a  set  of  five  storms.  His  fore¬ 
cast  results  for  these  five  storms  compared  frvorably  with  the  official  NHC  forecasts 
for  the  same  storms.  Curry  also  developed  theoretical  expressions  for  “confidence 
eUipses”  as  functions  of  the  model  coefficients  which  probabilistically  describe  the 
area  where  a  storm’s  future  position  is  likely  to  be.  These  are  discussed  in  following 
chapters. 

2.S  MoU’a  Model 

Following  Curry’s  threshold  ARMA  approach,  Mott  developed  a  trivariate 
model  that  added  maximum  sustained  wind  speed  as  an  explanatory  variable  (12). 
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In  each  latitude  band,  the  model  is: 


Xti4t  —  LAf—i  = 

4" 

+ 

+ 

LOt  —  LAt—i  = 
+ 
+ 
+ 

WSt-WSt-i  = 

+ 

+ 


—  LAt-i)  +  •  •  •  +  —  LAt-t) 

~  LOt-i)  +  •  •  •  +  —  LOt-t) 

MWSt-i  ~  WSt.i)  +  • . .  +  -  WSt-t) 

Cl  +  ttit 

^i,i{^At-i  —  LAt-2)  +  •  •  •  +  ^i,i{LAt-t  —  LAt-t) 

—  LOt-i)  4-  *  — I-  ^a,5(I'0<-s  —  LOf-e) 

MWSt.i  -  WS1.2)  4-  •  •  •  +  -  WSt.^) 

C2  4-  Ujt 

^i.iil^At^i  —  LAi-i)  4-  •  •  *  4-  ^i,t{LAt-i  —  LAt^n) 
<h2,i{LOt~i  —  LOt-i)  4-  •  •  •  4-  —  LOt-t) 

-  WS1.2)  4- .  •  •  4-  MWSi.2  -  WSi-t) 

C’a  4-  03« 


(7) 


(8) 


(9) 


where  LA  and  LO  are  latitude  and  longitude,  W5  is  wind  speed,  ^’s  are  the  model 
coefficients,  C’s  are  constants  and  the  a  tmms  are  random  errors. 

Starting  with  the  complete  data  set  through  1989,  Mott  reduced  it  similarly  to 
Curry  by  eliminating  storms  that  occurred  before  1945,  and  positions  observed  after 
landfall.  A  departure  from  Curry  was  retaining  data  of  tropical  depression  strength. 
This  was  useful  for  developing  a  model  for  all  stages  of  storms  since  one  cannot 
precisely  forecast  when  a  storm  will  change  above  or  below  tropical  storm  strength. 
It  also  allows  forecasts  to  be  made  earlier  in  a  storm’s  development.  Including  this 
data  provides  confidence  that  the  resulting  model  also  applies  to  both  strengthening 
and  weakening  storms. 
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i.S.t  Staiionarity.  Since  a  goal  of  Mott’s  research  was  to  create  a  moddi 
valid  for  all  hurricanes,  he  was  unwilling  to  limit  his  efforts  to  only  one  stationarity 
category  as  Curry  did.  Further,  Mott  noticed  most  storms  rarely  accderated  for 
very  long  and  ultimately  returned  to  what  Curry  termed  vdodty  stationary.  In  other 
words,  it  was  difficult  to  objectively  categorise  storms  into  one  of  Curry’s  stationarity 
categories  since  storms  appeared  to  change  category  over  time.  He  assumed  that  the 
probability  of  a  storm’s  fluctuation  from  its  mean  level  was  the  same  at  all  times, 
and  this  allowed  the  devdopment  of  a  model  with  fixed  coefficients  to  be  devdoped, 
regardless  of  stationarity  category. 

2.  S.  2  ModeUng  Scheme.  Mott  constructed  a  variety  of  univariate,  bivariate, 
and  tiivariate  models  on  which  he  based  his  coefficient  calculations.  He  used  SAS 
statistical  software  routines,  induding  backward  and  stepwise  regression,  to  calculate 
the  modd  coeffidents  and  their  statistical  significance.  Then,  each  modd  was  used 
to  forecast  storm  motion  for  each  data  point  and  great  drde  distance  error  statistics 
were  collected.  The  results  indicated  that  the  best  modd  was  trivariate,  using  past 
latitude,  longitude,  and  wind  speed  information  to  forecast  each  component.  Mott’s 
final  modd  and  coeffidents  are  also  induded  in  Appendix  B. 

2.3.3  OvUiera.  Both  Curry  and  Mott  noted  that  anomalous  tracks  occa¬ 
sionally  occurred  and  that  these  tracks  gave  the  greatest  forecast  errors.  Both  were 
concerned  about  the  possible  bias  this  data  introduced  into  their  respective  models 
because  the  behavior  of  these  storms  might  not  be  representative  of  typical  hurricane 
tracks.  Part  of  this  research  is  devoted  to  determining  identifying  characteristics  of 
these  storm  tracks  ac  ^  is  discussed  in  following  chapters. 
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2.4  Other  Modeb 


Other  approaches  to  hurricane  modding  and  forecasting  are  documented  in 
the  meteorological  literature.  Different  statistical  approaches  have  been  applied  to 
developing  hurricane  forecasting  models,  but  the  majority  of  models  use  meteoro¬ 
logical  input,  some  on  a  very  large  scale.  This  section  briefly  discusses  some  of  these 
other  modeb  and  their  relevance  to  a  time  series  approach. 

S.4.i  SiatuUcal  Modeb.  In  “An  Analysis  of  the  Error  Characterutics  of 
Atlantic  Tropical  Cyclone  Track  Prediction  Modeb”,  James  Kroll  examines  a  number 
of  hurricane  forecast  modeb  and  their  characterutics  (7).  He  develops  a  Combined 
Confidence  Weighted  Forecast  model  that  assigns  relative  weights  to  a  number  of 
existing  forecast  modeb  based  on  their  forecast  errors  (7:43-44).  Thu  amounts  to  the 
development  of  a  linear  program  for  minimiaing  the  aggregate  error  of  the  (NEC) 
modeb.  He  also  attempts  to  develop  linear  regression  modeb  using  a  number  of 
explanatory  or  predictor  variables,  such  as  vortidty,  storm  speed,  and  Julian  date, 
among  others  to  forecast  hurricane  position.  He  concludes  that  these  factors  are 
not  linearly  related  to  storm  track  (7:99);  that  is,  a  linear  regression  model  using 
these  factors  does  not  produce  good  track  forecasts.  Thu  reinforces  our  rationale  for 
developing  a  time  series  model. 

Merrill  develops  a  linear  regression  model  to  predict  tropical  cyclone  inten¬ 
sity  changes  (11).  Currently,  NEC  uses  SHIFOR  (iStatistical  Hurricane  Aitensity 
Forecast),  a  statistical  model  that  incorporates  initial  position,  date,  location,  mo¬ 
tion,  and  intensity  change  of  a  storm  to  predict  future  intensity  changes.  He  tests 
variables  such  as  sea  surface  temperature,  wind  information  at  various  leveb,  Julian 
date,  and  previous  intensity  changes  for  significance.  He  also  tests  “custom”  prod¬ 
ucts  of  the  variables  to  see  if  they  indicate  whether  non-linear  predictors  of  intensity 
change  are  useful.  Merrill’s  results  offer  slight  improvement  over  the  current  SHI¬ 
FOR  intensity  model,  but  they  are  not  statistically  significant  (11:24).  He  concludes 
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tropical  cyclone  intensity  changes  are  influenced  by  environmental  conditions,  but 
the  rdationships  are  weak  and  their  nsefolness  as  objective  forecasting  aids  is  limited 
(11:23).  This  shows  that  linear  regression  modeb  are  also  not  useful  for  forecasting 
change  in  tropical  cyclone  intensity. 

The  STRIKPA  model  explained  by  Jarrell  generates  a  hurricane  forecast  bul¬ 
letin  using  three  difficulty  classes  (5).  Instantaneous  and  time  integrated  land¬ 
fall  probabilities  are  generated,  analogous  to  the  STRIKPA  model  used  for  Pacific 
storms.  Use  of  these  difficulty  classes  is  unclear  and  the  reference  explaining  them 
is  unavailable,  so  our  ability  to  assess  the  utility  of  this  model  is  hampered. 

In  his  thesis,  “A  Kalman  Filter  With  Smoothing  for  Hurricane  Tracking  and 
Prediction”,  Asim  Mutaf  presents  a  Kalman  filter  model  to  forecast  smoothed  hur¬ 
ricane  tracks  (13).  This  method  was  also  applied  to  hurricane  wind  speed,  but 
necessary  wind  speed  data  was  frequently  unavailable  with  the  position  data  which 
resulted  in  a  frequently  unstable  model  (13:46).  Model  results  were  favorable  for 
the  two  control  storms  presented,  but  a  larger  data  base  is  necessary  to  validate  this 
technique.  Mutaf  mentions  this  method  is  particularly  useful  for  storms  with  scarce 
previous  track  data,  as  might  happen  in  some  Pacific  storms. 

2.4-2  Meteorohgical  Models.  The  vast  majority  of  hurricane  prediction 
modeb  use  large  amounts  of  observed  and  modeled  weather  data  and  meteorological 
equations.  The  NHC  generates  and  tests  most  of  these  for  possible  use  in  operational 
hurricane  forecasting.  Currently,  they  use  seven  different  modeb  simultaneously  to 
forecast  a  storm,  then  a  human  forecaster  subjectively  blends  these  forecasts  to  arrive 
at  the  official  forecast.  This  is  explained  in  “Modeb  for  the  Prediction  of  Tropical 
Cyclone  Motion  over  the  North  Atlantic:  An  Operational  Evaluation”  by  Neumann 
and  Pelissier  (17).  Since  this  report,  a  new  model,  NHC90,  has  been  developed. 
It  incorporates  a  variety  of  meteorological  information  including  climatological  and 
persistence  data,  and  observed  and  forecast  deep  layer  geopotential  heights  (similar 
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to  the  heights  of  piessure  levels).  It  is  based  on  a  previous  model  devebped  in 
1983  called  NHC83.  Daring  the  trial  phase  of  NHC83,  a  number  of  problems  in  the 
model  were  identified.  These  were  corrected  and  incorporated  into  the  new  model 
NHC90.  Input,  output  and  various  performance  considerations  are  described  in  "A 
Revised  National  Hurricane  Center  NHC83  Model  (NHCfiO)”  (16).  Unfortunately, 
the  researchers  devised  a  new  way  of  classifying  storms  into  north  and  south  sones 
that  obscures  similarities  with  previous  results.  As  well,  the  period  of  the  test  data 
is  differmit  from  other  research. 

FVanklin  and  DeMaria  indicate  the  importance  of  using  accurate  and  extensive 
wind  and  pressure  data  for  hurricane  forecasting  in  “The  Impact  of  Omega  Drop- 
windsonde  Observations  on  Barotropic  Hurricane  Forecasts”  (4).  Dropwindsondes 
consist  of  mostly  the  same  instruments  as  found  on  weather  balloons  (without  the 
balloon)  that  are  dropped  from  aircraft  near  or  in  hurricanes.  During  their  descent 
to  the  earth,  they  transmit  real-time  meteorological  observations  that  can  be  input 
into  forecast  models.  Barotropic  models  are  based  on  a  simplified  atmospheric  struc¬ 
ture  where  temperature  and  pressure  gradients  are  colinear.  Some  barotropic  models 
are  used  in  the  NHG  modeling  scheme.  More  extensive  wind  and  pressure  input  for 
these  models,  available  from  dropwindsondes,  markedly  improves  their  performance. 
A  significant  forecast  improvement  (10  to  15  %  reduction  in  error)  resulted  from  the 
use  of  dxopwindsonde  data,  but  the  authors  note  that  there  is  currently  no  theoreti¬ 
cal  basis  to  optimise  the  data  collection  pattern.  They  also  make  the  point  that  the 
cost  of  getting  the  data  is  more  than  offset  by  increased  forecast  accuracy  and  resul¬ 
tant  decreased  impact  possible  from  fewer  false  alarms  and  unnecessary  evacuations 
(4:390). 

Kurihara,  et  al,  describe  the  application  of  a  movable  mesh  model  to  forecast 
the  1985  hurricane  Gloria  (8).  Mesh  refers  to  an  arbitrarily  placed  grid  pattern  that 
specifies  points  on  the  earth’s  surface  where  weather  information  is  interpolated.  A 
movable  mesh  uses  the  same  grid  at  progressively  different  locations,  based  on  storm 
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movement.  A  fine  mesh  with  a  grid  spacing  of  and  two  coarse  mesha  of  I"  and 
l|‘’  were  used.  A  combination  of  coarse  and  fine  mesh  models  was  used  to  capture 
and  forecast  elements  of  the  storm  with  favorable  results  (8:2189).  However,  an 
important  caveat  is  brought  up,  namely,  the  selection  and  application  of  which  mesh 
and  when  to  use  it  is  highly  subjective  and  has  a  large  impact  on  the  results.  As 
well,  this  approach  has  been  applied  only  to  a  single  storm  after  the  fact  and  can’t 
yet  be  considered  to  be  generally  valid  for  other  storms. 

£.5  Compariaon  of  Pubhahed  Results 

Curry,  Mott,  and  NHC  use  great  circle  distance  to  report  errors  between  ob¬ 
served  and  forecast  positions.  It  can  be  computed  from: 

d  =  tQcoa~^[ain{LAf)ain{LA^  4-  coa{LAf)coa{LAt^coa{LOo  —  LOf)\  (10) 

where  d  is  in  nautical  miles,  subscript  /  denotes  the  forecast  obtained  from  the  model 
and  subscript  o  denotes  the  observed  position  from  the  best  track  data.  The  results 
are  shown  in  Teible  1  in  nautical  miles  and  corresponding  kilometers  in  parentheses. 

NHC  and  Curry  used  data  from  storms  of  tropical  storm  strength  or  greater. 
Mott  also  included  data  from  the  tropical  depression  stage.  Note  that  the  periods 
of  analysis  are  also  different.  As  a  restdt,  a  different  combination  of  storms  was 
used  by  each  source.  Curry  and  NHC  divide  the  number  of  initial  storm  positions 
approximately  in  half  by  creating  north  and  south  zones  such  that  roughly  as  many 
storms  start  above  the  dividing  latitude  as  below.  NHC  also  reports  model  error 
trends  using  a  rotated  coordinate  system.  Errors  are  further  characterized  by  NHC 
SIS  slow  or  fast,  and  right  or  left  of  track  using  track-oriented  coordinates  (17). 
Comparable  summaries  were  not  available  for  use  with  the  Curry  or  Mott  models 
and  so  are  not  reported  here. 
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At  first  ^ance,  it  would  appear  from  Table  1  that  the  forecasts  generated  by 
the  NHC  ate  the  best  since  their  average  reported  errors  over  all  regions  are  generally 
amaller  than  those  reported  by  Curry  or  Mott.  For  the  same  reasons,  Curry’s  model 
appears  to  perform  slightly  better  than  Mott’s.  Unfortunately,  the  differences  in  the 
data  sets  used  in  computing  the  these  averages  make  such  conclusions  tenuous  at 
best.  In  particular,  NHC  and  Curry  use  data  from  storms  of  tropical  strength  or 
greater,  while  Mott  includes  tropical  depresnons  as  wdl.  Perhaps  more  significantly, 
the  periods  of  analysis  are  also  different,  that  is,  different  combinations  of  storms 
are  used  by  each  source. 

The  effect  of  including  tropical  depressions  seems  to  inflate  the  average  errors. 
This  is  shown  by  Mott,  who  also  computes  forecast  error  for  Curry’s  model  using 
his  own  (unfiltered)  data  base.  As  seen  in  Table  2,  average  overall  forecast  errors 
for  Curry’s  model  are  larger  when  it  is  applied  to  all  storm  stages  from  1945  to  1983 
(including  tropical  depressions  as  computed  by  Mott)  than  when  tropical  depressions 
are  excluded  (as  computed  by  Curry).  The  implication  seems  to  be  that  tropical 
depressions  are  more  unpredictable  than  their  stronger  counterparts.  We  study  this 
in  more  detail  in  the  next  chapter. 


Forecast  Period  (Hours) 

12 

24 

48 

72 

Computed  by  Curry 

109.8 

214.6 

312.0 

Computed  by  Mott 

48.7 

269.8 

406.7 

Table  2.  Average  Overall  Forecast  Errors  for  Curry’s  Model 


Mott’s  results  should  also,  perhaps,  be  viewed  somewhat  suspiciously  as  evi¬ 
denced  in  Table  3,  which  summarizes  Mott’s  published  results  over  a  ntunber  of  time 
periods.  Clearly,  something  is  amiss  since  some  of  the  reported  averages  from  1945 
to  1988  are  smaller  than  either  of  the  corresponding  averages  computed  from  1945 
to  1983  and  1984  to  1988.  One  would  expect  the  1945  to  1988  average  to  be  between 
the  others. 
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Table  3.  Average  Overall  Forecast  Errors  with  IVopical  Depressions  (computed  by 
Mott) 
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III.  Hurricane  Prediction  Regions 


This  chapter  describes  the  methods  and  rationale  by  which  we  attempt  to  fur¬ 
ther  Mott’s  time  series  hurricane  forecast  model  to  include  probabilistic  prediction 
regions  for  future  storm  position.  Since  Mott’s  final  model  is  the  basis  of  all  manip¬ 
ulations  and  calculations,  a  first  step  is  to  compare  the  difference  in  forecast  errors 
between  using  tropical  depression  data  in  the  model,  as  Mott  did,  and  removing 
this  data,  as  NHC  does.  Once  the  effect  of  this  is  established,  the  transition  from  a 
simple,  rectangular  prediction  region  to  a  bivariate  elliptical  region  is  developed.  A 
Monte  Carlo  simulation  is  used  to  generate  these  probability  regions  and  the  utility 
of  this  method  is  demonstrated. 

S.l  Tropical  Depression  Data 

One  of  Mott’s  goals  was  to  develop  a  model  that  was  useful  for  forecasting 
aU  stages  of  tropical  cyclones,  regardless  of  strength  or  stationarity  category.  He 
noted  that  it  was  nearly  impossible  to  tell  when  or  for  how  long  a  particular  storm 
would  (or  wouldn’t)  increase  or  decrease  in  intensity.  Thus,  he  favored  including 
all  data  from  the  storms  in  the  database,  regardless  of  wind  speed.  Mott  included 
data  from  tropical  depressions,  that  is,  storms  and  storm  segments  vrith  sustained 
winds  of  less  than  34  knots  in  developing  his  models.  Because  they  are  generally 
less  organized,  their  motion  can  be  expected  to  be  influenced  more  by  external 
meteorological  conditions  than  are  stronger  tropical  storms  and  hurricanes  (2:314). 

On  the  other  hsmd,  no  other  researchers  include  tropical  depression  data  in 
their  model  development  or  forecast  error  reports,  and  we  vrished  to  check  if  including 
this  data  may  have  biased  both  Mott’s  model  and  his  reported  error  statistics.  This 
data  would  need  to  be  removed  and  the  model  coefficients  recalculated  if  Mott’s 
results  were  statistically  different  without  tropical  depression  data. 


Position  reports  from  storms  with  wind  speeds  less  than  34  knots  thus  were 
filtered  from  the  storm  database.  All  data  greater  than  tropical  storm  strength  was 
retained;  only  points  with  wind  speed  below  34  knots  were  removed.  Mott’s  final 
model  was  then  run  first  with  all  the  original  data,  then  vrith  the  filtered  data.  Error 
statistics  were  generated  for  each  database  and  compared.  Results  from  the  model 
with  the  unfiltered  data  are  shown  in  Table  4  and  results  with  the  filtered  data  are 
shown  in  Table  5.  Differences  are  slight,  and  hypothesis  tests  confirm  the  differences 
are  statistically  insignificant  at  the  99%  level. 

These  results  imply  that  the  moddl  coefficients  apply  equally  well  to  tropical 
depressions,  tropical  storms,  or  hurricanes  and  that  inclusion  of  the  tropical  depres¬ 
sion  data  does  not  substantially  bias  Mott’s  reported  measures  of  forecast  accuracy. 
Since  we  agree  with  Mott  that  this  data  ought  to  be  included  in  developing  a  viable 
model,  this  research  proceeded  without  recomputing  the  model  coefficients. 

S.2  Prediction  Region  Development 

An  advantage  that  statistical  models  have  over  meteorological  models  is  that 
they  naturally  provide  a  means  of  describing  a  region  within  which  the  model  would 
predict  a  storm  would  likely  be.  Assuming  that  the  random  errors  contained  in 
our  model  have  a  bivariate  normal  distribution,  it  follows  that  the  errors  associated 
with  a  storm’s  forecast  position  also  have  a  bivariate  normal  distribution  and  the 
latitude  and  longitude  of  its  future  position  are  each  normally  distributed.  One  can 
then  easily  construct  sepuate  lOO(l-a)  %  prediction  intervals  for  a  storm’s  future 
latitude  and  longitude  if  the  variances  of  the  forecewt  errors  are  known  or  can  be 
estimated.  In  particular,  if  we  let  LAt^i  and  LOt+i  denote  the  forecasted  latitude 
and  longitude  at  I  time  steps  in  the  future,  and  let  (r\  and  <7,  denote  the  variances 
of  the  /-step  ahead  forecast  errors  in  the  latitude  and  longitude  directions,  then 
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OVERALL  FORECAST  ERROR  SUMMARY 


FORECAST 

#OBS 

MEAN 

STDEV 

WS  MEAN 

WS  STDEV 

6HR 

7069. 

13.8837 

22.1548 

0.1944 

5.8290 

6734. 

41.4018 

39.2187 

0.2803 

9.8451 

24HR 

6063. 

103.5616 

78.0337 

-0.0430 

15.5701 

48HR 

4743. 

243.2267 

163.4678 

-1.9407 

21.7397 

72HR 

3639. 

376.3317 

239.5433 

-3.9957 

FORECAST  ERRORS  FOR  LAT  BAND  1 


FORECAST 

#OBS 

MEAN 

STDEV 

WS  MEAN 

WS  STDEV 

6HR 

573. 

6.6948 

12.2494 

0.1686 

5.1662 

12HR 

560. 

23.7970 

22.0288 

0.2158 

9.1464 

24HR 

539. 

62.5049 

52.7147 

0.1717 

48HR 

495. 

152.6815 

108.ie97 

-1.3243 

72HR 

454. 

251.3975 

165.0699 

-3.2103 

30.3096 

FORECAST  ERRORS  FOR  LAT  BAND  2  j 

FORECAST 

#OBS 

MEAN 

STDEV 

WS  MEAN 

WS  STDEV 

6HR 

1093. 

9.7225 

26.8267 

0.3743 

6.7017 

12HR 

1064. 

31.3281 

39.3789 

0.7524 

11.6103 

24HR 

1008. 

77.7159 

56.5139 

1.6058 

18.4564 

48HR 

909. 

186.4410 

103.7929 

1.4624 

24.9959 

72HR 

820. 

302.0040 

166.7822 

-0.0922 

27.5448 

FORECAST  ERRORS  FOR  LAT  BAND  3 

FORECAST 

#OBS 

MEAN 

STDEV 

WS  MEAN 

WS  STDEV 

6HR 

1325. 

11.5459 

18.9466 

-0.1292 

6.1360 

12HR 

1302. 

35.0921 

35.6572 

-0.3335 

10.8593 

24HR 

1253. 

87.7139 

69.4099 

-1.1214 

17.6346 

48HR 

1114. 

mm 

154.8390 

-2.8742 

23.5915 

72HR 

925. 

366.1059 

250.1919 

-3.8596 

25.7184 

FORECAST  ERRORS  FOR  LAT  BAND  4 

FORECAST 

#OBS 

STDEV 

WS  MEAN 

WS  STDEV 

6HR 

1485. 

12.6875 

18.8105 

0.3545 

6.8699 

12HR 

1423. 

39.9202 

0.5239 

11.0885 

24HR 

1280. 

106.5221 

74.0204 

-0.2781 

16.1690 

48HR 

1028. 

270.3296 

168.5778 

-3.8992 

20.4066 

72HR 

784. 

426.6084 

242.8052 

-7.7132 

20.3057 
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1  FORECAST  ERRORS  FOR  LAT  BAND  5 

FORECAST 

#OBS 

MEAN 

STDEV 

WS  MEAN 

WS  STDEV 

6HR 

1212. 

14.8686 

19.9300 

0.1229 

4.5284 

12HR 

1176. 

47.7177 

37.2576 

0.1059 

7.4863 

24HR 

1108. 

127.2018 

81.5491 

11.5795 

48HR 

825. 

311.9072 

174.6172 

-2.0081 

15.0116 

72HR 

515. 

499.6970 

252.0175 

-4.4721 

17.3789 

FORECAST  ERRORS  FOR  LAT  BAND  6 

FORECAST 

#OBS 

MEAN 

STDEV 

WS  MEAN 

WS  STDEV 

6HR 

946. 

20.4669 

25.1238 

0.2415 

4.4064 

12HR 

900. 

58.7448 

47.8595 

7.2818 

24HR 

705. 

141.7330 

94.0775 

0.0713 

11.1526 

48HR 

334. 

333.3453 

185.2412 

-3.1191 

15.9542 

72HR 

129. 

524.4794 

271.5908 

-9.1660 

18.1830 

FORECAST  ERRORS  FOR  LAT  BAND  7 

FORECAST 

#OBS 

MEAN 

STDEV 

WS  MEAN 

WS  STDEV 

6HR 

435. 

27.9528 

27.7191 

0.3133 

5.3766 

12HR 

309. 

66.8525 

43.9893 

0.2426 

8.0306 

24HR 

170. 

169.1263 

87.1928 

-0.6208 

11.6383 

48HR 

38. 

445.2230 

245.4738 

0.7943 

12.8953 

72HR 

12. 

798.6640 

306.7503 

7.9552 

10.9006 

Table  4.  Forecast  Errors  With  Tropical  Depression  Data 
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FORECAST 

#OBS 

6HR 

1471. 

OVERALL  FORECAST  ERROR 


MEAN 


40.1329 


100.4248 


233.9071 


380.6521 


24HR 


48HR 


72HR 


1385. 


1209. 


877. 


654. 


SUMMARY 


STDEV  WS  MEAN 


19.0288  0.2165 


36.3064  0.0662 


73.5527  -0.9955 


151.9349  -5.6711 


231.5491  -8.9159 


WS  STDEV 


6.4126 


10.0711 


19.4610 


20.9687 


FORECAST  ERRORS  FOR  LAT  BAND  1 


FORECAST 

#OBS 

MEAN 

STDEV 

WS  MEAN 

WS  STDEV 

6HR 

5.4433 

11.3793 

-0.2193 

3.7339 

12HR 

106. 

22.9768 

23.1243 

-0.6172 

6.1661 

24HR 

103. 

64.2772 

48.7063 

-2.0682 

9.6871 

48HR 

97. 

178.2417 

106.7414 

-6.5209 

14.9612 

72HR 

94. 

313.8691 

155.1456 

-10.2317 

17.0372 

FORECAST  ERRORS  FOR  LAT  BAND  2 


FORECAST 

#OBS 

MEAN 

STDEV 

6HR 

191. 

8.5425 

13.8736 

12HR 

180. 

29.9391 

23.6592 

24HR 

162. 

76.4328 

51.2605 

48HR 

130. 

189.4273 

102.6489 

72HR 

115. 

294.7719 

130.8438 

:mean 


.5217 


-1.3900 


-11.1588 


-18.9672 


WS  STDEV 


7.3472 


12.1232 


19.2433 


22.0947 


FORECAST  ERRORS  FOR  LAT  BAND  3 


FORECAST 

#OBS 

MEAN 

STDEV 

WS  MEAN 

WS  STDEV 

6HR 

258. 

10.7047 

15.5602 

-0.8213 

4.9106 

12HR 

252. 

34.4744 

28.1402 

-1.9214 

9.5103 

24HR 

237. 

87.4906 

55.9169 

-4.0593 

18.5607 

48HR 

195. 

fmm 

137.2368 

-9.2227 

22.9932 

72HR 

159. 

372.8253 

231.7108 

-8.6047 

23.5791 

FORECAST  ERRORS  FOR  LAT  BAND 


FORECAST 

#OBS 

MEAN 

STDEV 

WS  MEAN 

WS  STDEV 

6HR 

335. 

10.5417 

16.0678 

0.3530 

12HR 

321. 

36.8703 

31.6514 

0.3739 

24HR 

279. 

97.5051 

69.1440 

-0.9616 

18.1356 

48HR 

225. 

249.4264 

171.8977 

-3.6698 

19.5612 

72HR 

169. 

432.9576 

264.6495 

-3.7649 

19.2547 
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FORECAST  ERRORS  FOR  LAT  BAND  5 


FORECAST 

#OBS 

MEAN 

STDEV 

WS  MEAN 

WS  STDEV 

6HR 

276. 

12.7738 

19.3964 

0.6836 

4.6113 

12HR 

262. 

44.6530 

37.0867 

1.0515 

6.5987 

24HR 

246. 

119.4761 

79.2823 

1.5729 

9.9131 

48HR 

176. 

281.8597 

m 

14.1480 

72HR 

105. 

448.7641 

267.2599 

17.0769 

1  FORECAST  ERRORS  FOR  LAT  BAND  6  ] 

FORECAST 

#OBS 

MEAN 

STDEV 

WS  MEAN 

WS  STDEV 

6HR 

200. 

18.0333 

25.6591 

0.3508 

4.4162 

12HR 

192. 

54.6750 

48.5064 

0.6810 

7.5500 

24HR 

149. 

133.1868 

97.1556 

0.5706 

11.5582 

48HR 

51. 

267.9384 

167.6344 

•4.9871 

14.8084 

72HR 

12. 

497.8853 

261.2467 

-19.7256 

20.9371 

FORECAST  ERRORS  FOR  LAT  BAND  7 

FORECAST 

#OBS 

MEAN 

STDEV 

WS  MEAN 

WS  STDEV 

6HR 

104. 

25.3356 

25.1876 

0.7421 

5.4312 

12HR 

72. 

69.9979 

48.7447 

0.9876 

7.4625 

24HR 

33. 

158.6578 

85.7521 

4.7570 

48HR 

3. 

343.6791 

11.8737 

2.3850 

3.4518 

Table  5.  Forecast  Errors  Without  Tropical  Depression  Data 
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lOO(l-a)  %  piedictioa  intervals  on  a  storm’s  future  position  will  be  given  by 


LAt+l  ±  Za/ttTl 
LOt^i  ±  Sa/a<ra 

where  Za/a  is  a  critical  value  from  a  standard  normal  distribution  corresponding  to 
the  specified  significance  level,  a.  If  latitude  and  longitude  were  independent,  one 
would  then  have  a  rectangular  prediction  region  with  probability  of  100(1  —  a)’%  of 
a  forecast  falling  inside  of  it. 

Unfortunately,  latitude  and  longitude  can’t  be  considered  independent  for  this 
purpose  because  future  storm  coordinates  are  constrained  by  the  physical  process  of 
storm  motion.  Clearly,  for  any  given  latitude,  a  storm  can’t  randomly  move  to  any 
longitude  and  vice  versa.  Additional  effort  is  needed  to  calculate  a  prediction  region. 
In  particular,  according  to  standard  statistical  results  [see,  for  example,  Johnson  and 
Wichem  (6)],  if  a  storm’s  future  latitude  and  longptude  at  (-position  reports  ahead 
have  a  bivariate  normal  distribution,  we  can  write: 

1  I  ^  ^IrOt-n  —  LOt+i'^ 

—  LAt+i\  ( LOt+i  -  LOt+i 

<ri  J  \  <r% 

where  LAt^i  and  LOt+i  zxe  the  storm’s  future  latitude  and  longitude,  p  is  the  corre¬ 
lation  between  them,  <ri  and  0*3  are  their  respective  standard  deviations,  LAt+i  and 
LOt+i  are  their  respective  means,  and  is  the  critical  value  from  a  Chi-square 
distribution  with  two  degrees  of  freedom  at  the  desired  significance  level. 

3.2.1  Curry’s  V*  Method.  Curry  develops  a  series  of  equations  to  define  a 
prediction  region  based  on  his  so-called  ^  form  of  the  ARMA  model  equations  in 


) 


<xi«  =l-a 


(11) 
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hit  diuertatioii  (3).  In  order  to  calculate  the  prediction  ellipset,  the  model  equa¬ 
tions  need  to  be  manipulated  into  this  form.  This  involves  algebraically  rewriting 
his  bivariate  AR  models  as  infinite-order  MA  models  whose  coefficients  are  denoted 
by  ^*8.  Once  these  are  derived  for  the  model  equations  within  each  latitude  band, 
expressions  for  the  desired  variances  and  correlation  can  be  developed.  Within  each 
latitude  band,  these  computations  are  rather  tedious,  but  straightforward.  They 
become  much  more  complicated  when  the  storm’s  forecasted  track  crosses  latitude 
bands  and  require  a  complicated  updating  scheme  to  be  used.  In  particular,  when 
forecast  positions  cross  a  threshold,  ijf  values  of  the  previous  band  are  used  for  obser¬ 
vations  and  forecasts  in  the  old  bsmd  and  values  of  the  new  ^r’s  are  used  for  forecasts 
in  the  new  band,  and  the  computations  become  somewhat  intractable  (3:79). 

Curry  uses  the  resultant  expressions  to  generate  confidence  ellipse  axes  (3:68- 
75).  The  equations  are  based  on  the  one-step  ahead  forecast  errors.  These  errors  are 
assumed  to  have  a  multivariate  normal  distribution  because  each  is  random  ‘Vhite 
noise”  error.  Curry’s  resulting  prediction  ellipse  has  essentially  the  same  form  as 
equation  (11),  but  erroneously  uses  a  critical  value  from  a  normal  distribution,  rather 
than  the  appropriate  Chi-square  distribution  (6:126). 

In  order  to  apply  these  results  to  Mott’s  model,  an  arduous  reworking  of  the 
equations  would  be  required  to  account  for  three  explanatory  variables,  instead  of 
the  two,  latitude  and  longitude,  that  Curry  uses.  This  would  be  time  consuming  to 
develop  and  check,  u  well  as  code.  We  thus  consider  a  more  appealing  and  intuitive 
approach. 

S.8.2  Monte  Carlo  Simulation.  An  alternative  to  the  deterministic  com¬ 
putation  of  the  desired  variances  and  correlation  coefficient  is  to  run  the  model 
numerous  times  with  simulated  random  errors  to  generate  the  statistical  equivalent 
to  the  solution  of  Curry’s  ^  equations.  This  is  conceptually  simpler,  quicker  to  code, 
and  urives  at  approximately  the  same  results.  Sample  statistics  are  used  to  esti- 
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mate  tlie  needed  variances  and  correlation  coefficient.  It  requires  relatively  greater 
amounts  of  processing  to  obtain  solutions,  but  is  still  a  reasonable  approach  to  use 
operationally. 

This  method  proceeds  by  generating  a  number  of  simulated  storm  tracks  using 
the  model’s  defining  equations  (given  by  Equations  7-9)  with  simulated  values  from  a 
normal  distribution  used  to  represent  the  random  model  errors.  Within  each  latitude 
band,  these  model  errors  are  each  assumed  to  have  a  mean  of  sero  and  a  variance 
equal  to  the  mean  squared  error  found  by  Mott  in  fitting  his  model  equations.  The 
correlation  between  the  latitude  and  longitude  errors  is  assumed  to  be  sero  since 
Curry’s  results  suggest  this  is  a  reasonable  first  approximation  (3).  (Curry’s  results 
show  that  the  correlations  within  each  latitude  band  are  small.  In  most  cases,  we 
would  fail  to  reject  the  hypothesis  that  the  correlation  was  sero  at,  say,  a  5  % 
level  of  significance.  Because  the  correlations  are  all  relatively  small,  it  is  doubtful 
that  accounting  for  them  explicitly  would  substantially  alter  the  resulting  prediction 
regions.)  The  FORTRAN  code  to  accomplish  this  is  included  in  Appendix  C. 

In  order  to  produce  a  prediction  ellipse,  the  1-period  ahead  position  from  each 
simulated  storm  track  is  recorded  and  these  are  then  used  to  estimate  the  variances  of 
the  /-position  ahead  latitude  and  longitude  forecast  errors  as  well  as  the  correlation 
between  them.  These,  in  turn,  are  used  in  defining  the  equation  of  the  corresponding 
ellipse. 

In  our  code,  200  simulated  storm  tracks  are  generated  each  time  a  prediction 
ellipse  is  desired.  If  the  number  of  such  tracks  is  large  (as  we  assume  in  this  case), 
then  the  sample  variances  will  approximate  the  actual  variances  and  the  use  of 
Chi-square  critical  values  remains  valid.  Otherwise,  criticzd  values  from  Hotelling’s- 
t  statistic  should  be  used  (6).  Figure  3  shows  an  example  of  the  24-hour  99% 
forecast  ellipse  and  corresponding  probability  rectangle  for  Hurricane  Hugo  (1989) 
approximately  24  hours  before  landfall. 
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S.S.S  Landfaii  Prediction  Regions.  The  idea  of  landfall  prediction  is  to 
determine  the  probability  that  a  storm  will  make  landfall  in  a  given  section  of  coast. 
An  obvious  approach  based  on  our  Monte  Carlo  simulation  is  to  run  the  model 
numerous  times  from  a  position  offshore  and  keep  track  of  the  proportion  of  tracks 
that  cross  a  given  section  of  coast.  This  proportion  provides  an  estimate  of  the 
probability  that  a  storm  will  make  landfall  in  that  region.  Generally,  this  will  enable 
us  to  assess  the  location  of  land&ll,  but  not  its  exact  timing,  however.  Predicting  the 
time  of  landfall  (within  a  margin  of  a  few  hours)  is  not  as  important  as  predicting 
the  area  affected  for  effective  evacuation. 

Another  approach  is  to  note  the  coastal  area  "covered”  when  a  prediction 
ellipse  intersects  the  coast.  Varying  the  signiiicance  level,  a,  will  generate  contours 
of  equal  probability  which  intersect  different  sections  of  coast  and  provide  an  idea  of 
the  likelihood  of  landfall  at  a  particular  location. 

The  more  accurate  the  a  model  is,  the  more  accurate  the  forecast  it  gives.  The 
next  chapter  presents  our  efforts  to  improve  the  time  series  model  by  attempting  to 
identify  potential  outliers  in  the  data. 
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IV.  Examination  of  Anomalous  Storms 

This  chapter  presents  the  results  of  attempts  made  to  identify  possible  mod¬ 
ifications  to  Mott’s  model  that  could  potentially  be  made  to  improve  its  forecast 
accuracy. 

4.1  Identifying  Anomaloiu  Siorma 

Curry  and  Mott  commented  on  the  possible  impact  of  anomalous,  or  outlier, 
data  on  the  statistical  results  of  their  respective  models.  Outliers  introduce  poten¬ 
tially  significant  hidden  bias  into  the  modds  because  it  is  difficult  to  determine  if 
a  particular  storm  will  generally  fit  the  model  ahead  of  time.  This  research  effort 
is  aimed  at  identifying  and  characterizing  apparently  anomalous  storms  in  order  to 
determine  if  they  have  any  common  characteristics  which  could  be  exploited  to  im¬ 
prove  the  forecast  model.  The  storms  were  ranked  according  to  forecast  error,  then 
resulting  patterns  were  analyzed. 

4.1.1  Ranking  Forecast  Errors.  One  approach  to  identifying  outliers  is 
to  rank  the  forecasts  by  magnitude  of  error  and  identify  storms  with  the  largest 
errors.  For  example,  storms  with  forecast  errors  of  greater  than,  say,  three  standard 
deviations  above  the  mean  forecast  error  could  be  classified  as  potentially  anomalous 
since  one  normally  expects  only  about  1%  of  all  forecast  errors  to  exceed  this  value. 
If  outliers  could  be  identified  early  in  their  life,  they  could  be  modeled  separately. 
Our  intent  is  to  eventually  remove  outliers  from  the  data  base  and  recompute  the 
model  coefficients.  The  resulting  model  would  potentially  have  greater  accuracy  by 
removing  the  bias  introduced  by  anomalous  storms. 

In  this  research,  storms  were  ranked  according  to  the  magnitude  of  their  av¬ 
erage  and  maximum  24-  and  72-hour  forecast  errors.  After  ranking,  the  tracks  of 
some  storms  with  an  error  greater  than  three  standard  deviations  above  the  mean 
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were  plotted  to  view  their  behavior.  Moet  outliers  exhibited  eccentric  tracks,  some 
looping  and  others  accelerating  beyond  the  model’s  ability  to  predict  them.  Others 
tended  to  move,  often  rapidly,  from  west  to  east  while  the  storms  that  were  forecast 
best  tended  to  move  steadily.  Some  of  these  tracks  are  plotted  in  Figure  2. 

Tables  6  and  7  show  the  list  of  storms  ranked  by  magnitude  of  forecast  error. 
The  format  shows  the  ID  number  of  the  storm  from  the  data  base,  the  month  and 
year  of  occurrence,  maximum  error  during  the  storm,  and  average  error  during  the 
storm.  A  cursory  inspection,  in  response  to  Mott’s  conjecture  that  outliers  may 
tend  to  be  either  early  season  or  late  season  storms,  shows  most  storms  with  poor 
forecasts  occurred  after  August.  This,  however,  is  not  consistent  since  quite  a  few 
good  forecasts  occur  after  August,  also.  As  well,  the  peak  of  storm  activity  occurs 
in  September,  so  it  should  not  surprising  that  most  of  the  worst  forecast  errors  are 
experienced  during  that  time  (18:203). 

Histograms  of  the  average  and  maximum  errors  were  prepared  in  an  effort  to 
identify  if  there  were  any  potential  outliers  that  could  be  better  modeled  separately. 
These  are  shown  as  Figures  3  through  6.  There  did  not  seem  to  be  a  clear  trend 
with  the  data  available.  In  fact,  the  errors  seemed  to  approximately  follow  normal 
distributions.  The  expected  pattern  produced  by  outliers  would  show  up  as  a  peak 
further  from  the  main  concentration  of  errors,  indicating  a  distinct  group  of  forecast 
errors.  The  results  show  no  such  patterns,  even  when  the  resolution  (range  of  errors) 
is  changed.  Another  distinguishing  characteristic  might  be  those  forecasts  that  lie 
outside  of  the  prediction  ellipse  developed  in  the  next  section.  That  is,  if  a  storm’s 
porition  were  to  fall  outside  the  99%  probability  region,  it  could  be  considered  a 
potential  outlier  and  similar  analysis  could  be  done.  Time  constraints  did  not  allow 
this  to  be  accomplished  for  this  thesis. 

Other  patterns  of  errors  are  not  intuitively  forthcoming.  Evidently,  factors 
other  than  those  measured  in  the  database  are  likely  involved.  It  seems  reasonable 
to  assume  that  some  meteorological  factor  or  process  drives  outlier  behavior  that 
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cannot  be  modeled  with  data  from  the  current  database.  With  this  result,  our 
research  ended. 
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Storm  ID 

Storm  Date 

Max  Erjor 

Avg  Error 

775 

Oct  1975 

406.37 

123.81 

847 

Oct  1947 

160.28 

169 

Aug  1969 

144.95 

651 

Sep  1951 

] 

158.36 

856 

Nov  1956 

4  11.24 

230.84 

581 

Sep  1981 

411.46 

172.21 

276 

Aug  1976 

412.49 

150.08 

572 

Sep  1972 

414.51 

162.61 

752 

Oct  1952 

417.06 

114.14 

350 

Aug  1950 

422.50 

115.05 

668 

Sep  1968 

430.35 

215.19 

854 

Oct  1954 

431.25 

115.28 

561 

Sep  1961 

431.64 

98.38 

354 

Aug  1954 

436.80 

103.74 

959 

Oct  1959 

443.21 

155.97 

461 

Sep  1961 

459.01 

100.50 

258 

Aug  1958 

462.57 

152.64 

851 

Oct  1951 

469.30 

122.44 

1061 

Nov  1961 

474.20 

209.11 

980 

Oct  1980 

479.23 

157.52 

849 

Sep  1949 

490.48 

143.88 

871 

Sep  1971 

501.67 

85.97 

1278 

Nov  1978 

515.69 

252.87 

470 

Aug  1970 

532.77 

250.77 

955 

Sep  1955 

546.17 

116.77 

1250 

Oct  1950 

604.81 

157.60 

481 

Aug  1981 

1205.31 

182.48 

Table  6.  24  Hour  Ranked  Errors 
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Storm  ID 

Storm  Date 

Max  Error 

Avg  Error 

778 

Sep  1978 

1016.37 

451.82 

1253 

Oct  1953 

1024.41 

737.40 

258 

Aug  1958 

1030.62 

530.74 

1058 

Oct  1958 

1032.94 

415.17 

367 

Sep  1967 

1043.50 

377.25 

450 

Sep  1950 

1045.85 

319.91 

958 

Sep  1958 

1050.38 

645.49 

553 

Sep  1953 

1060.15 

691.55 

572 

Sep  1972 

1061.29 

817.39 

651 

Sep  1951 

1130.37 

580.13 

172 

May  1972 

1136.55 

812.74 

650 

S?p  1950 

1188.04 

377.25 

1166 

Nov  1966 

1188.35 

790.57 

1061 

Nov  1961 

1200.92 

717.42 

467 

Sep  1967 

1206.40 

677.32 

665 

Oct  1965 

1226.99 

514.65 

673 

Sep  1973 

1227.87 

568.29 

149 

Aug  1949 

1240.87 

949.23 

748 

Sep  1948 

1282.29 

578.41 

561 

Sep  1961 

1285.91 

403.93 

863 

Oct  1963 

1288.61 

525.64 

1278 

Nov  1978 

1296.03 

838.33 

849 

Sep  1949 

1347.05 

538.78 

1145 

Oct  1945 

1347.68 

1005.65 

1469 

Oct  1969 

1369.30 

605.66 

653 

Sep  1953 

1423.60 

1235.82 

350 

Sep  1950 

1461.37 

452.54 

856 

Nov  1956 

1471.26 

970.60 

871 

Sep  1971 

1486.30 

470.33 

1255 

Oct  1955 

1501.15 

1064.36 

980 

Oct  1980 

1507.28 

_ 

410.28 

Table  7.  72  Hour  Ranked  Errors 
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Figure  2.  Anomalous  Storm  Tracks 
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Figure  3.  Histogram  of  24  Hour  Average  Forecast  Error 
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Figure  4.  Histogram  of  24  Hour  Maximum  Forecast  Error 
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Figure  5.  Histogram  of  72  Hour  Average  Forecast  Error 
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Figure  6.  Histogram  of  72  Hour  Maximum  Forecast  Error 


V.  Summary  and  Conclusions 

Cuiy  initiated  research  on  using  a  time  series  model  to  describe  hurricane 
motion  (3).  Forecasts  made  with  his  model  compare  favorably  with  forecasts  made 
by  the  NHC.  Curry  also  devdoped  a  method  of  generating  prediction  ellipses  as  a 
means  of  describing  an  area  where  a  storm  might  move,  rather  than  just  a  potentially 
misleading  point  forecast. 

Mott  explored  and  improved  the  model;  the  best  forecasts  were  made  with  the 
addition  of  a  third  parameter,  wind  speed. 

A  number  of  other  researchers  have  added  to  the  field  of  hurricane  modding  and 
prediction.  Two  sources  using  statistical  approaches  point  out  that  linear  regression 
models  are  not  effective  for  predicting  hurricane  motion  or  intensity,  even  with  a 
variety  of  variables  (7)  (11).  Most  other  researchers  have  explored  meteorological 
aspects  to  varying  degrees.  Despite  better  computers  and  more  extensive  input 
data,  the  result  of  meteorological  research  seems  to  be  a  gradually  decreasing  rate 
of  improvement  in  forecasting  ability,  as  measured  by  forecast  errors. 

Both  Curry  and  Mott  suggest  that  their  modek  could  be  biased  by  outlier 
data  and  that  identification  of  potential  outliers  could  lead  to  forecast  improvement, 
possibly  through  the  devdopment  of  a  separate  modd  to  handle  outlier  storms. 
Unfortunatdy,  while  outliers  can  be  easily  identified  after  the  fact,  there  are  as  yet 
no  objective  criteria  to  apply  as  a  storm  occurs.  Despite  our  efforts,  we  were  unable 
to  improve  the  situation,  although  we  did  establish  that  time  of  season  does  not  seem 
to  be  a  consistent  indicator.  This  is  an  area  where  meteorological  experts  could  apply 
thdr  experience  to  suggest  likdy  criteria  to  explore.  Being  able  to  separate  outliers 
from  “normal”  storms  should  increase  the  accuracy  of  the  time  series  modd. 

A  time  series  modd  is  ultimatdy  limited  because  it  does  not  account  for  the 
underlying  meteorological  processes  driving  hurricane  devdopment.  But  the  surpris¬ 
ing  accuracy  with  such  a  simplistic  modd  suggests  it  could  have  even  better  residts 
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when  oonpied  with  n  metecwologicdi  bans.  Althoo^  there  are  opportunities  for  fine 
tuning  the  time  series  performance,  snch  as  increasing  the  number  of  latitude  bands, 
the  gmtest  improvement  may  lie  in  exploring  meteordogical  avenues. 

Devdopnmnt  of  prediction  ellipses,  especially  for  landfall  prediction  was  one  of 
our  goab.  Curry’s  ^  method  would  provide  a  straightforward  solution,  except  some 
of  the  quantities  needed,  namely  covariances,  are  not  easily  obtained  or  estimated. 
Applying  the  ^  method  to  Mott’s  model  would  be  arduous  and  prone  to  errors.  The 
Monte  Carlo  approach  was  conceptually  easier  to  accomplish  and  the  results  are 
materially  the  same  as  from  Curry’s  method.  (One  change  from  Curry’s  formulation 
is  our  use  of  a  Chi-square  critical  value,  rather  than  one  from  a  standard  normal 
distribution.)  Using  Mott’s  final  modd  for  greatest  applicability,  prediction  ellipses 
give  a  quantitative  measure  of  probability  for  a  position  forecast. 

This  approach  can  be  used  to  generate  landfall  prediction  zones  when  com¬ 
bined  with  fine  resolution  coastUne  information.  The  result  is  greater  confidence  in 
the  bounds  of  an  area  at  risk  and  attendant  opportunity  to  minimize  unnecessary 
evacuations. 
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Appendix  A.  Determining  Tropical  Cyclone  Position  and  Intensity 

Prom  Satellite  Imagery 

This  appendix  describes  procedures  used  to  determine  tropical  cyclone  posi¬ 
tion  and  intensity  firom  satellite  imagery.  This  material  is  based  on  NO  A  A  Technical 
Report  NESDIS  11  (14)  and  a  National  Weather  Association  Meteorological  Mono¬ 
graph  (15). 

Tropical  meteorologists  have  been  using  satellite  images  for  monitoring  tropical 
storms  for  almost  25  years.  Satellite  imag^  is  important,  especially  at  early  storm 
development  stages  and  in  remote  areas  where  weather  measurements  are  otherwise 
unavailable.  Satellites  continue  to  carry  more  and  better  instruments  as  time  goes 
on.  Infrared  (IR)  and  visible  light  sensors  provide  the  primary  sources  of  imagery, 
but  microwave  images  and  image  processing  will  shortly  provide  enhanced  capability 
and  a  greater  variety  of  information. 

Visual  imagery  is  used  whenever  possible  because  of  the  relatively  good  resolu¬ 
tion  and  imaging  characteristics.  Infrared  images  are  important  because  they  enable 
the  analyst  to  monitor  a  storm’s  activity  continuously  night  and  day.  The  IR  image 
of  a  tropical  cyclone  appears  similar  to  that  seen  in  a  visual  image  of  the  same  storm. 
In  fact,  the  day  by  day  changes  in  the  cloud  pattern  can  often  be  followed  in  the 
same  way  they  are  in  visual  imagery. 

To  determine  reliable  center  coordinates  from  visible  and  IR  imagery  on  a 
routine  basis,  a  systematic  approach  is  necessary.  Such  an  approach  is  outlined  in 
the  foUowing  sections.  It  consists  of  five  steps  which  are  designed  to  help  the  anal3rst 
avoid  common  errors  in  tropical  cyclone  positioning  when  using  satellite  imagery. 

The  five  steps  are: 

1.  Locate  the  overall  pattern  center 

2.  Examine  for  small  scale  features 
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3.  compaie  center  location  with  forecast 

4.  Compare  center  with  previous  pattern  center 

5.  Make  final  center  adjustments 

The  first  step  is  to  determine  the  center  of  the  overall  cloud  pattern  of  the 
tropical  cyclone.  This  “first  guess”  estimate  of  the  center  is  determined  by  one  of 
two  methods:  either  by  locating  the  focal  point  of  the  broad  scale  cloud  line  and 
band  curvature  of  the  system,  or  by  comparing  the  cloud  pattern  being  analysed 
with  patterns  in  the  model  of  tropical  cyclone  development  which  includes  the  center 
placement  for  each  pattern. 

In  the  second  step  of  the  analysis  procedure,  the  area  around  the  pattern  center 
found  in  step  1  is  examined  for  small  scale  cloud  feature  indications  that  often  give 
an  improved  center  location.  The  third  step  is  to  compare  the  center  arrived  at  in 
the  first  two  steps  with  the  extrapolated  track  (or  forecast)  position  for  the  time  of 
the  satellite  image.  This  step  brings  to  Ught  the  possible  errors  due  to  wrong  cloud 
center  selection  or  due  to  displacements  caused  by  vertical  wind  shear.  It  presumes 
that  short  period  forecast  positions  are  fttirly  accurate. 

In  the  fourth  step,  the  center  placement  in  a  previous  storm  and  its  relation  to 
cloud  features  is  compared  with  the  current  center  placement  relative  to  similar  cloud 
features.  This  step  is  important  when  step  three  reveals  a  possible  location  error. 
The  fifth  step  involves  possible  adjustments  to  the  center  location  when  indicated 
by  either  gridding  errors,  by  high  satellite  viewing  angles,  or  by  the  displacement  of 
the  cloud  system  center  away  from  the  sur&ce  wind  center  of  the  cyclone. 

The  confidence  one  places  in  the  center  location  depends  on  four  primary  fac- 
tors.  They  are:  the  type  of  cloud  features  defining  the  center,  the  vertical  integrity 
of  the  vortex,  the  grid  accuracy,  and  the  type  of  imagery  used  in  the  analysis.  When 
using  enlarged  visible  pictures  with  accurate  gridding,  the  following  guidelines  are 
recommended: 


1)  For  low  level  cloud  lines  defining  a  tight  center  and  for  well-defined  eyes, 
the  accuracy  is  comparable  to  that  of  reconnaissance  locations  (±0.2°  latitude). 

2)  For  tightly  curved  upper  level  features  and  moderately  defined  low  level 
features,  accuracies  of  ±0.3°  latitude  are  expected. 

3)  For  multiple  center  patterns  and  sheared  patterns,  average  errors  of  0.4° 
latitude  or  24  nmi  are  to  be  expected.  For  this  last  category  of  patterns,  it  should 
be  remembered  that  multiple  wind  and  pressure  centers  may  also  involve  errors  in 
the  verification  data. 

4)  The  confidence  ratings  listed  above  in  2)  and  3)  should  be  lowered  for  IR 
(without  eyes),  and  for  cloud  systems  beginning  to  show  the  effects  of  strong  vertical 
shear. 

Additional  technical  details  are  contained  in  the  references  and  they  should  be 
consulted  for  further  information. 
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Appendix  B.  The  Time  Series  Models 


B.l  Chary  *8  Model 

Curry  found  the  second,  third,  and  fifth  coefficients  were  small,  so  he  dropped 
them  from  his  model.  The  following  tables  show  the  coefficients  he  reported,  sepa:- 
rated  by  latitude  band. 
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B.S  Mott’s  Model 


Equations  7-9  show  the  form  of  Mott’s  model.  The  following  table  gives  the 
^  coefBdents  by  latitude  band.  The  first  row  lists  the  five  latitude-latitude  (^i,n) 
coefficients,  the  second  row  lists  the  latitude-longitude  coefficients,  the  third 

row  lists  the  latitude-wind  speed  (^3,n)  coefficients.  The  fourth  through  sixth  rows 
list  the  longitude—  •  •  coefficients  and  the  seventh  through  ninth  rows  list  the  wind 
speed—  •  •  coefficients. 
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Appendix  C.  Monte  Carlo  FORTRAN  Code  to  Generate  Probability 

Ellipse 

C  NOITE.FOR 

c 

C  A  PROGRAM  TO  GENERATE  HURRICAIE  PROBABILITY 
C  ELLIPSES  USING  MONTE  CARLO  SIMULATION 
C 

INTEGER  J,  K.  I,  M.  NSTEP,  NPOINT,  L 

REAL  LATDEV.LONDEV.Al ,A2,A3.LAL0RH0,C0VSUN.DISC 

REAL  C0EF(16,14) .C0EF2(ie,7) ,C0V13(7) ,C0V23(7) 

REAL  LATSUN,LONSUM,HSSUN.LATRND,LONRND,VSRND.CHI 
REAL  LATMEAN,HSNEAN.LONMEAN,LATVAR,LONVAR,LASTEP 
REAL  RH012(7),  RH013(7),  RH023(7) 

REAL  F0RE(1000,60),  PLOTLAT(IOOO) >  PLOTLON(IOOO) 

REAL  LAL0C0V,C0V12(7) ,LATNSE(7) ,L0NMSE(7) ,NSNSE(7) 

REALMS  SEED 

CHARACTER*20  OLDFIL,  NEHFL 

C  *****  MEAN  SQUARED  ERRORS  FROM  MOTT’S  RESULTS 

DATA  LATMSE/ . 025 , . 033 , . 079, . 066 , . 087 , . 124 , . 239/ 

DATA  LONMSE/ . 046 , . 071 , . 206 , . 112 , . 137 , . 192 , . 479/ 

DATA  WSMSE/21. 43, 36. 21 ,40. 67,34. 53, 18. 41, 18. 7, 23. 43/ 

DATA  C0V12/7*0./ 

DATA  C0V13/7*O./ 

DATA  C0V23/7*0./ 

SEED»1234567 

MP0IMT*6 
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c  *****  OPEI  FILE  FOR  OUTPUT  DATA 
HRITE(*,250) 

250  FORIIAT(X,*IIPUT  FILEIANE  FOR  ELLIPSE  DATA’) 

READ(*,’(A20)O  lEHFL 

0PEI(3O,  FILE-IEHFL,  STATUS-'HEVO 

C  *****  6ETTIIG  6  DATA  POIHTS  TO  FORECAST  FROM 

DO  680  I>0,6 

«RITE(*,670)  I 

READ  ♦,  F0RE(l,7-iy 

WIITE(*,671)  I 

READ  *,  F0RE(1,27-I) 

raiTE(*,672)  I 

READ  *,  F0RE(1,47-I) 

WRITE(30,739)  I 

WRITE(30,740)  F0RE(1 ,7-1) ,F0RE(1, 27-1) ,F0RE(1, 47-1) 
DO  675  J>2,HP0IHT 
FORE(J ,7-I)=F0RE(l ,7-1) 

F0RE(J,27-I)=F0RE(1,27-I) 

F0RE(J,47-I)=F0RE(1 ,47-1) 

675  COHTIHUE 
680  COHTIHUE 
HRITE(30,741) 

670  F0R1!AT(X,'IHPUT  LATITUDE  AT  TIME  T-*,!!) 

671  FORMAT(X,»IHPUT  LOHGITUDE  AT  TIME  T-»,I1) 

672  FORMAT(X,»IHPUT  WIHD  SPEED  AT  TIME  T-',I1) 

739  FORMAT(X,' OBSERVED  IHPUTS  FOR  TIME  T-»,I1) 
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740  F0RIfAT(3(2X,F6.2)) 

741  FORMAK*  ») 

C  *****  GET  THE  MATRIX  THAT  HAS  THE  FORECASTIIG  COEFFICIEMTS 
HRITE  (*,501) 

501  FORMAT (X.* FILE  lAME  OF  MATRIX  VITH  THE  FORECASTIIG  COEFFICIEITS? ' ) 

READ(*.>(A20)*)  OLDFIL 

OPEI (UlIT-24 , FILE-OLDFIL . STATUS^  >  OLD ’ , lOSTAT-IERROR . ERR»500) 

00  510  M  -  1,16 

READ(24,*,EID«520)  C0EF(M,1),  C0EF(M,2),  C0EF(M,3), 

+  C0EF(M,4).  C0EF(M,5),  C0EF(M,6),  C0EF(M,7),  C0EF(M,8), 

C0EF(M,9),  C0EF(M,10),  C0EF(M,11),  C0EF(M,12),  C0EF(M,13), 

+  C0EF(M,14) 

510  CONTIMUE 
520  CLOSE (24) 

GO  TO  900 

500  HRITE(*,5O7)IERR0R 

607  FORMAT (»♦♦*  CAIHOT  OPEI  FILE  ♦♦♦’,18) 

c  *****  get  coefficiehts  for  PREDICTIHG  HS 

900  HRITE  (*,901) 

901  F0RMAT(1X,’FILE  lAME  OF  MATRIX  WITH  THE  HS  COEFFICIEHTS?’) 

REAO(*,  ’(A20)’)  OLDFIL 

OPEN (UIIT»44 , FILE»OLDFIL , STATUS= ’ OLD ’ , IOSTAT=IERROR , ERR=999) 

DO  910  M  »  1,16 

READ(44,*,END»920)  C0EF2(M,1),  C0EF2(H,2),  C0EF2(M,3), 

+  C0EF2(M,4),  C0EF2(M,6),  C0EF2(M,6),  C0EF2(M,7) 

910  COHTINUE 
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920  CLOSE (44) 
GO  TO  922 


999  WRITE(*,907)IERROR 

907  FORMAT CAIIOT  OPEl  HS  FILE  ****,18) 

C  *****  GENERATE  SIMULATED  VALUES 
C  *****  (HSTEP  SIX  HOUR  PERIODS  AHEAD) 

922  LATSUM*>0 

L0ISUM»O 

HSSUM«0 

RH012 (K) -C0V12 (K) / (LATMSE(K) *LOIMSE(K) ) 

RH013 (K) »C0V13 (K) / (LATMSE(K) *HSMSE(K) ) 

RH023 (K) »C0V23 (K) / (LOHMSE (K) *HSMSE(K) ) 

ISTEP-4 

DO  300  I>1,NP0IHT 
DO  560  J-1,ISTEP 

C  *****  get  PROPER  LATITUDE  BAND  FOR  FORECAST  MODEL 
K»0 

IF(F0RE(I,6+J).LT.16.0)K=1 

IF (  (FOREd , 6+J)  . LT .  20 . 0)  .  AID .  (FOREd , 6+ J)  .  GE .  15 . 0)  )K*2 
IF ( (FORECI ,6+J) .LT . 25 . 0) . AID . (FORECI ,6+J) . GE . 20 . 0) )K=3 
IF ( (FORE (I , 6+J ) . LT . 30 . 0) . AID . (FORE (I , 6+ J) . GE . 25 . 0) )K=4 
IF ( (FORE ( I , 6+ J ) . LT . 35 . 0) . AID . (FORE (I , 6+ J) . GE . 30 . 0) ) K»5 
IF ( (FORE ( I , 6+ J) . LT . 40 . 0) . AID . (FORE (I , 6+ J) . GE . 35 . 0) )K=6 
IF(F0RE(I,6+J) .GE. 40.0)1=7 
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C  *****  RilDOH  ERROR  GEIERATIOI 
R-RAIDOM(SEED) 

S-RAIDOH(SEED) 

T-RODOMCSEED) 

LATRMD- (R** . 1349- (1-R) ** . 1349) / . 1975 
LOIRMD- (S** . 1349- (1-S) ** . 1349) / . 1975 
WSRMD- (T** . 1349- (1-T) ** . 1349) / . 1975 
A1«LATRID*LATIISE(K) 

A2-LATRlD*L0IMSE(K)*RH012(K)-fLQIR8D*L0IIfSE(K)*SQRT(l-RHQ12(K)**2) 
A3-LATRID*WSIfSE(K) *RHQ13(K) +LQHRn}*VSMSE(K) * (RH023 (K) -RH012 (K) 

+  *RHQ13(K))/SQRT(l-RH012(K)**2)-i-VSRHD*HSMSE(K)*SQRT((l 

+  -RH012 (K) **2) * (1-RH013 (K) **2) - (RH023 (K) -RH012 (K) 

+  *RH013 (K) ) **2) /SQRT(1-RH012 (K) **2) 

C  ««***  forecast  latitude  USIIG  liPUTTED  FORECAST  MODEL  COEFFICIEITS 
FOREd  ,7+J)=«F0RE(1 ,6+J) 

++COEF  ( 1 ,  K)  *  (FORE  ( 1 , 6+ J) -FOREd ,  5+ J)  ) 

++C0EF(2  ,K)  *  (FOREd  ,5+J)-F0REd  ,4+J)  ) 
++C0EF(3,K)*(F0RE(I,4+J)-F0RE(I,3+J)) 

++COEF (4 , K) * (FORE (I , 3+ J) -FORE(I , 2+ J) ) 

++COEF (5 , K) * (FORE ( I , 2+ J) -FORE(I , 1+ J) ) 

++COEF  (6 ,  K)  *  (FORE  ( 1 , 26+ J)  -FOREd ,  25+ J )  ) 

++COEF  (7  ,K)  *  (FOREd ,  25+ J)  -FORE(1 , 24+J)  ) 

++C0EF(8 ,  K)  ♦  (FORE  (1 , 24+J)  -FOREd ,  23+ J)  ) 

++COEF  (9  ,K)  *  (FOREd ,  23+J)  -FORE(1 , 22+J)  ) 
++C0EF(lO,K)*(F0RE(I,22+J)-F0RE(I,21+J)) 

++COEF ( 1 1 ,K) ♦ (FORE (I , 46+J) ) 

++C0EF(12,K)*(F0RE(I,45+J)) 


53 


•MCOEFda  ,E)  *(FORE(1 ,44-t-J)  ) 

++C0EF(14,K)*(F0RE(I,43+J)) 

+<i-C0EF(16  .K)  *(FORE(1 .42-t-J)  ) 

•t^4COEF(16,K)+Al 

C  *****  FORECAST  LOIGITUDE  USIIG  IIPUTTED  FORECAST  MODEL  COEFFICIEITS 
FORECI . 27+ J) -FORE( I , 26+ J) 

++COEF  (1 ,  R+7)  *  (FOREd ,  6+J)  -FORECI ,  5+J)  ) 

++COEF (2 , K+7) * (FORECI , 6+J) -FORE(I ,4+J) ) 

++COEF (3 , K+7) * (FORECI , 4+J) -FORECI ,3+ J) ) 

++COEF (4 , K+7)* (FORECI .3+J) -FORECI ,2+J) ) 

++C0EF(6 ,I+7)*(F0RE(I , 2+J)-F0RE(I, 1+J)) 

++COEF (6 ,K+7) ♦ (FORECI , 26+J) -FORE(I , 25+J) ) 

++COEF (7 , K+7) * (FORECI , 26+ J) -FORE(1 , 24+J) ) 

++COEF  (8 , K+7)*  (FORECI ,  24+J) -FORECI ,  23+J)  ) 

++COEF (9 , K+7) * (FORECI , 23+ J) -FORECI , 22+ J) ) 
++C0EF(10,K+7)*(F0RECI,22+J)-F0RE(I,21+J)) 

++COEF(ll ,K+7) * (FORECI ,46+ J) ) 

++COEF ( 12 , K+7) * (FORECI , 46+ J) ) 

++C0EF(13 ,K+7) * (FORECI ,44+ J) ) 

++COEF ( 14 , K+7) * (FORECI , 43+ J) ) 

++C0EF(16 ,K+7) ♦ (FORECI ,42+ J) ) 

++C0EF(16,K+7)+A2 

C  *****  FORECAST  VIID  SPEED  USIIG  IIPUTTED  FORECAST  MODEL  COEFFICIEITS 
FORECI, 47+J)-COEF2(l,K)*(FORE(I,6+J)-FORE(I, 6+J)) 
++C0EF2(2,K)*(F0RE(I,6+J)-F0RE(I,4+J)) 

++C0EF2 (3 ,K) * (FORECI ,4+J) -FOHE(I ,3+J) ) 
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4-k:QEF2  (4 .1)  *(FORE(1 ,3+J) -FORECI .  2-i-J)  ) 

■f4COEF2  (5  .K)  *  (FORE  (1 . 2+ J) -FOREd .  UJ)  ) 

+4C0EF2  (6 ,  K)  *  (FOREd .  26-»^J)  -FORECI .  26+J)  ) 

++COEF2  (7  ,K)  •(FOREd  ,25+J)-F0REd ,  24+ J)  ) 

++C0EF2 (8 ,K) ♦ (FORE (I , 24+J) -FORE(I , 23+ J) ) 

++C0EF2  (9,1)*  (FOREd  .23+J)  -FOREd ,  22+J)  ) 

++C0EF2  (10 ,  E)  •  (FOREd .  22+J)  -FORE(1 , 21+ J)  ) 
++C0EF2  (11  .!)•  (FOREd  ,46+J)  ) 

++C0EF2  (12 .  K)  *  (FOREd .  45+ J)  ) 

++C0EF2  (13  .K)  *  (FOREd  .44+ J)  ) 

++C0EF2 (14 . K) *(FORE(I .43+ J) ) 

++C0EF2 ( 1 5 . K) * (FORE ( I . 42+ J) ) 

++C0EF2(16.I)+A3 

VRITE(30.729)  I.J 

WRITE(30.730)  F0RE(I.7+J)  .FOREd. 27+J)  .FOREd. 47+J) , 
+  A1.A2.A3 

miITE(30.741) 

560  corrilUE 
LATSUN-LATSUN+FORE (I . 7+ISTEP) 

LORSUM»LOISUN+FORE(I . 27+ISTEP) 

VSSUM-HSSUM+FORE (I . 47+ISTEP) 

HRITE(+.710)  LATSUN.LOHSUH.VSSUl! 

300  CORTIRUE 
miITE(30.741) 

710  F0RMAT(3(X.F8.2)) 

729  FORMATd.’AT  RPOIRT  ’.I2.»  ARD  RSTEP  ».I2. 

+  ^  FORECAST  VALUES  ARE:’) 
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730  F0RIUT(X,*UT-»,F6.2,‘  L0I-*,F6.2,*  »IID-',F6.2, 

*  RAIDOM  ERRORS  ARE\2(2X,F6.4).21,F7. 3) 

C  *****  CALCUUTII6  SIMPLE  DEVIATIOIS,  COVARIAICE, 

C  *****  AID  CORREUTIOI  OF  FORECAST  DATA 

LATNEAI-LATSUM/IPOirr 

LOIMEAI-LOISUM/IPOIIT 

HSNEAI-VSSUM/IPOIIT 

HRITE(«.667) 

667  FORMATCX.’FIIISHED  HITH  FORECASTIIGO 
DO  578  L-l.lPOirr 

LATVAR-  (LATMEAI-FORECL .  7-flSTEP)  )  *«2-i>LATVAR 
LOIVAR-(LOniEAI>FORE(L ,27+ISTEP) )**24L0IVAR 
COVSUM-(LATMEAI-FORE(L  .7-t-ISTEP)  ) 

*(L0IMBAl-F0RE(L,27-t-ISTEP))'K:0VSUM 

578  COITIIUE 

HRITE(30,750)  LATMEAI.LOIHEAI, LATVAR, LOIVAR 
HRITE(30,741) 

UTDEV-  (LATVAR/IPOIIT)  ** .  5 
LORDEV- (LORVAR/IPOIIT) ** . 5 
LALOCOV-COVSUM/IPOirr 
LALORHO-LALOCOV/  (UTDEV*LOIDEV) 

LATVARaLATVAR/IPOIMT 

LOIVAR-LOIVAR/IPOIIT 

750  F0RMAT(X,'LATMEAI»»,F6.2,»  LOMMEAM-’ ,F6.2, »  LATVAR* » , 
+  F6.3,'  LOIVAR*’ ,F6.3) 

C  *****  POIHTS  01  ELLIPSE  CALCULATIOM 
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CHI-9.21 


L-0 

USTEP-0. 

219  DISC«(l.-UL0RH0**2)*(CHI-US'rEP*«2) 
PLOTLATCD-LATMEAI+USTEP^LATDEV 

PLOTLOI(L)«LOIHEAl4LOIDEV*  (LALORHO*LASTEP-(-SQRT(DISC)  ) 
L-L+1 

USTEP-USTEP+.2 
IF(LiSTEP*«2.LE.CHl)  GOTO  219 

220  LASTEP-USTEP>.2 
IF(UlSTEP*«2.GT.CHI)  goto  221 
DISC»(1 .  "LAI.0BH0**2)*(CHI-LASTEP**2) 
PLQTLAT(L)-UTlfEAI-<-LASTEP4>LATDEV 

PLQTLQM  (L)  -LOlHEAI-t-LOIDEV*  (LALORHO*LASTEP--SqRT(DISC)  ) 

L-L+1 

GOTO  220 

221  LASTEP-LASTEP+.2 
IFCLASTEP.GT.O.)  GOTO  222 

DISC- ( 1 , -LAL0Ra0**2) ♦ (CHI-LASTEP**2) 

PLOTLAT (L) -LATHEAI+LASTEP+LATDEV 

PLOTLOM(L)-LQIIIEAK+LOn)EV*(LALaRaO*LASTEP+SqRT(DISC)) 

L-L+1 

G0T0221 

C  mm***  OUTPUT  ELLIPSE  POIHTS  TO  A  FILE 

222  HRITE(30.720) 

DO  215  I-0,L-1 

WHTE(30,270)  PLOTLAT(I)  ,PLOTLOH(I) 
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215  COITIIUE 
CL0SE(3O) 

720  F0RI!AT(2X.'LlT\4X,a0lO 
270  F0R1UT(Z.2(X.F6.2)} 

EID 

FUICTIOI  RAMDOM(SEED) 

REAL*8  PROD,  SEMI.  SEED 
PROD-16807. DO*SEED 
SENl-DNOD (PROD , 2 147483647 . DO) 
RA1D0M-SEMI*0 .4656613E-9 
SEED-SEMI 
RETURI 

EID 
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